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ABSTRACT 


In  order  to  study  effects  of  velocity  and  acceleration  of  a  flow 
to  the  pressure  un  an  obstacle,  a  small  sphere  with  a  pressure  hole 
was  placed  in  a  periodically  pulsating  jet.  Both  the  instantaneous 
pressure  on  the  sphere  and  the  instantaneous  velocity  of  the  flow 
field  when  the  sphere  was  absent  were  measured.  By  using  periodic 
sampling  and  averaging  tcchiaques,  only  periodic  or  deterministic 
component  of  the  signal  was  recovered  and  any  random  component 
caused  by  turbulence  suppressed.  The  measured  surface  pressure 
was  expressed  in  terms  of  this  measured  velocity  and  acceleration 
of  the  flow.  A  simple  inviscid  theory  was  developed  and  the  experimental 
results  were  compared  with  it. 
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1.  INTRODUCTION 


Measurement  of  the  instantaneous  value  of  the  static  pressure 


inside  a  turbulent  flow  is  still  an  unresolved  problem  because  the 


pressure  probe  placed  in  the  unsteady  flow  field  represents  a  solid 


boundary  and  the  surface  pressure  at  a  point  always  contains  important 


contributions  from  the  inertial  effects  in  the  flow  around  the  body.  For 


steady  flows  it  is  relatively  easy  to  design  "static"  probes  that  read 


the  static  pressure  by  using  potential  theory  to  calculate  the  pressure 


distribution  around  the  body  of  the  probe,  but  the  same  locations  do 


not  give  instantaneous  static  pressure  in  imsteady  flow. 


In  order  to  understand  batter  the  behavior  of  small  probes  placed 


in  unsteady  flows,  an  experiment  was  performed  in  a  relatively  simple 


configuration.  A  small  sphere  was  placed  in  a  pulsating  jet  and  the 


instantaneous  surface  pressure  fluctuations  were  measured.  Since 


in  such  a  flow  there  is  a  strong  random  component  of  fluctuations  due 


to  the  turbulence,  a  special  signal  processing  technique,  namely  periodic 


sampling  and  averaging  was  performed  on  all  signals  in  order  to  enhance 


the  periodic  (or  deterministic)  component  and  to  suppress  the  random 


component. 


Furthermore,  in  order  to  provide  a  guide  for  the  assessment  of 


the  results  a  simple  inviscid  theory  was  developed  and  the  experinenlal 


results  were  compared  with  it. 
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2.  EXPERIMENTAL  FACILITY  AND  PROCEDTJRE 


The  time-dependent  flow  for  measurements,  the  pulsating  flow, 
is  produced  by  the  scheme  indicated  in  Fig-  1-  The  steady  air  flow 
provided  by  the  centrifugal  blower  is  divided  into  two  streams,  the 
one  flowing  through  the  duct  Dj  and  the  other  through  the  duct  D^- 
One  flow  component  is  discharged  through  the  duct  Dj  and  is  periodically 
iiitercepted  by  a  rotating  disc  with  16  holes.  As  a  result  periodic 
pressure  fluctuations  are  induced  also  in  the  other  branch,  D^,  through 
which  the  other  flow  component  is  discharged  and  pulsating  air  jet  is 
produced.  The  frequency  and  emplitude  of  the  velocity  fluctuation  in 
the  jet  can  be  varied  within  limits  varying  the  speed  of  the  rotating 
disc  and  the  stand-off  between  the  disc  and  the  nozzle  attached  to  Dj. 

The  pulsating  velocity  of  the  jet  was  measured  by  a  constant 
temperature  hot-wire  anemometer  followed  by  a  linearize r. 

Since  hot-wire  signal  as  obtained  from  the  outpiit  of  the  linearize r 
still  contains  random  fluctuations  due  to  turbulence  in  the  jet,  special 
signal  processing  technique,  periodic  sampling  and  averaging,  was 
performed  by  using  "Wavefrom  Eductor"  Model  TDH9  of  Princeton 
Applied  Research.  The  output  of  that  instrument  gives  the  ensemble 
average  of  a  large  number  (typically  1000  -  2000)  of  sweeps  and  the 
random  compor^nt  is  correspondingly  reduced  by  a  factor  of  V  1000 
V^OGO  or  30-40.  As  the  result  most  of  the  random  components  is 
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suppressed  and  only  periodic  or  deterministic  component  is  recovered. 

The  necessary  synchronizing  pulses  were  obtained  from  the  rotating 
disc  by  using  a  photo-cell  pick  up.  Fig.  2  shows  an  example  of  the 
periodic  velocity  fluctuation  measured  at  a  distance  x  =  6  cm  from 
the  nozzle  at  a  pulsating  frequency  of  f  =  450  Hz.  The  same  technique 
was  used  to  obtain  the  periodic  pressure  fluctuation  on  the  sphere  as 
described  below.  Schematic  diagram  of  the  instrumentation  is  shown  in 
Fig.  3. 

Three  spheres  of  the  same  diameter  but  with  a  pressure  hole  at 
different  locations  along  the  meridian  were  used  as  models.  They  are 
designated  as  S,  O  and  R  as  shown  in  Fig.  4.  The  diameter  of  the 
sphere  is  6.  35  mm  and  the  diameter  of  the  pressure  hole  is  0.  5  mm, 
located  at  the  meridian  angle  0  =  0°  (opposite  stem),  45°  and  90°  for 
S,  O,  and  R  respectively.  The  pressure  hole  is  connected  to  a  frequency 
modulated  condenser  micnphone  through  the  stem  as  shown  in  Fig.  4, 
and  the  pressure  fluctuations  is  converted  to  the  electrical  signal.  In 
earlier  period  of  the  experiment  a  transistiorized  version  of  the  circuit 
used  by  Einstein  and  Li  was  used  as  the  FM  oscillator  and  the  detector. 
In  later  period  an  improved  circuit  with  the  carrier  frequency  within  the 
range  of  commercial  FM  broadcast  was  built  into  ihe  probe  as  shown  in 
Fig.  4  and  an  FM  receiver  EICO  CORTINA  3200  was  used  as  the  detector. 

The  sphere  was  rotated  around  its  center  by  a  lathe  turn  labl»-  up 
to  a  maximum  angle  of  +  20°  (Fig.  5)  so  by  using  all  three  models  it 


was  possible  to  orient  the  pressure  holes  on  the  sphere  at  every  5 
between-20°  to  110°.  The  overall  static  sensitivity  was  found  to  be 
57  and  2.  4  ySl/y^oOT  respectively-  The  overall  frequency  characteristic-* 
of  the  pressure  probe  depends  on  the  diameter  and  length  of  the  pipe 
leading  to  the  microphone  as  well  as  on  the  volume  of  the  cavity  in 


front  of  the  microphone  diaphragm.  Fig-  6  shows  the  calculated 
frequency  characteristics  for  the  combination  of  a  pipe  of  diameter 
1.  0  mm  and  length  20.  5  mm  and  for  two  values  of  the  cavity  of  volume 

O 

1.  65  and  89.  1  mm  ,  respectively-  Resonances  occur  around  800  Hz 
and  2500  Hz  respectively.  The  dynamic  sensitivity  of  the  probe  would 
have  been  best  obtained  by  the  Coupler  Method,  which  is  based  on  a 
comparison  with  a  standard  microphone.  Unfortunately,  however,  no 
standard  microphone  was  available,  so  an  alternative  procedure  was  used 
by  utilizing  a  theoretical  relation  obtained  for  a  sphere.  According  to  the 
theoretical  calculation  given  in  a  later  section,  the  pressure  fluctuation 
on  a  sphere  is  made  up  of  two  terms,  one  being  proportional  to  the 
fluctuating  velocity  itself  and  the  other  one  proportional  to  the  time  derivative 
of  the  tluctuating  velocity.  At  the  forward  stagnation  point  the  non- 
dimensional  coefficients  of  the  two  terms  calculated  by  the  theory  arc  1.  0 
and  1.  5  respectively.  By  assuming  a  phase  lag  in  the  measured  pressure 
fluctuation  due  to  the  lead  pipe  and  the  cavity  volume,  the  ratio  of  the  two 
coefficientSjdynamic  sensitivityjfor  the  stagnation  point  was  calculated 
as  described  in  the  following  chapter  from  the  fluctuating  velocity  and  its 
time  derivative  obtained  by  a  hot-wire  placed  at  the  location  of  the  sphere 


center  but  in  the  absence  of  the  sphere,  TJ.o  phase  lag  so  determined 


varied  between  I  and  30°,  and  the  dynamics  sensitivity  was  found  to  be 
10%  -  40%  higher  than  the  static  sensitivity. 


3.  EXPERIMENTAL  RESULTS 

First  measurements  were  made  within  the  jet  which  is  discharged 
from  the  duct  D2  through  a  contracting  nozzle  to  the  ambient  air  (Fig.  1  )• 
for  a  mean  velocity  of  about  20  m/s  at  the  jet  axis.  Fig.  7  shows  the 
velocity  distribution  across  the  steady  non  pulsating  jet  at  two  stations, 

X  =  3  and  6  cm  dovr*Mrtrearn  from  the  nozzle.  The  core  of  the  jet, 
defined  as  the  region  in  which  the  velocity  is  greater  than  90%  of  the 
velocity  on  the  axis  (^o  ®  diameter  1.  55  cm  at  x  =  3  cm  and  1.28  cm 

at  6  cm,  respectively.  Fig.  8  shows  the  instantaneous  velocity  profiles 
across  the  fluctuating  jet  at  a  station  x  =  5  cm  for  a  frequency  of  450  Hz. 
The  core  diameter  defined  as  above  is  about  1. 20  cm,  which  is  a  little 
less  than  that  of  the  steady  jet.  Fig-  9  shows  the  streamwise  variation 
cf  the  amplitude  and  phase  velocity  of  the  velocity  fluctuation  along  the 
jet  axis.  Circles  are  described  to  indicate  the  location  and  size  of  the 
sphere.  When  the  sphere  is  placed  at  x  =  3  cm,  the  ratio  of  the  core 
diameter  of  the  steady  jet  to  the  sphere  diameter  is  2.  4,  but  the  amplitude 
of  velocity  fluctuation  still  varies  by  about  30  per  cent  over  the  region 
occupied  by  the  sphere.  When  the  sphere  is  placed  at  x  =  6  cm,  the 
ratio  of  the  diameters  is  reduced  to  2.  0,  but  the  variation  of  the  amplitude 
cf  the  velocity  fluctuation  becomes  negligible.  The  phase  velocity  determined 
by  two-hot  wires  separated  by  a  fixed  streamwise  distance  is  approximately 
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constant  over  the  region  covered  by  experiment.  The  actual  value  was 
13-2  m/s,  about  two  thirds  of  the  mean  velocity  at  the  jet  axis.  The 
phase  velocity  is  regarded  as  the  velocity  of  the  travelling  vortex  ring 
produced  at  the  nozzle  traveling  downstream. 

Fig.  10  shows  the  mean  (D.  C- )  pre  sure  distribution  on  the  sphere 
placed  at  X  =  3  cm  in  the  steady  jet-  The  calculated  value  are  also  shown. 

The  Reynolds  number  based  on  the  sphere  diaumeter  and  the  velocity  on  the 

(4  A  o 

axis  is  -j^-8.  8  x  10  -  As  shown  later  in  the  theory,  the  instantaneous 
surface  pressure  at  a  point  can  be  expressed  as 

P  =  P.  +  a|-u^  + 

where  U  is  the  instantaneous  streamwise  velocity  component^  a  is 
the  diameter  of  the  sphere,  A  and  B  are  the  coefficients  of 
the  contribution  from  the  instanteaneous  velocity  and  from  the 
instantaneous  acceleration  respectively. 

The  coefficients  A  and  5  are  strong  functions  of  the  coordinates 
(the  azimuth  angle)  and  they  may  also  weakly  depend  on  the  Reynolds 
number. 

For  the  periodically  pulsating  flow,  the  above  equation  can  be  written 
more  specifically  as 

p'fT)  =  A 

where  P^=(P  “  PoX/ ?  ,  T='t/'T  ,  T  =  Period 

for  a  particular  location  on  the  sphere.  After  the  periodic  sampling  and 
averaging,  values  of  all  oi  p(X)  ,  U^(T)  U('€)  were 


anc 


sampled  at  ten  equally  spaced  points  over  the  entire  period,  and  the 


coefficient  A  and  B  were  determined  by  the  method  of  least  square 


so  that  the  mean  square  error 


L  [  P(Ti)  -  tBQ-Ul-ri)jj 


becomes  minimum. 


Fig.  11  to  14  show  the  experimental  values  of  A  and  B,  the  tv  o  coefficients 


of  influence.  The  values  determined  on  the  three  models  S,  O  and  R  are 


indicated  in  these  figures  by  the  symbols  of  square,  circle  and  triangle. 


respectively.  Solid  curves  give  the  theoretical  prediction  based  on  the 


assumption  of  inviscid  irrotational  flow  described  in  the  following  chapter. 


4.  THEORY 


With  a  view  to  gaining  an  insight  into  the  problem,  a  theory  is 


developed  to  predict  the  pressure  acting  on  a  sphere  placed  in  a  time- 


dependent  flow  of  an  inviscid  fluid.  A  solution  of  the  problem  relevant 


to  the  subject  was  treated  by*  Lamb  as  follows.  The  pressure  p  on  a 


sphere  moving  with  time -dependent  velocity  V  in  an  infinite  mass  of 


fluid  at  rest  at  infinity  is  given  as 


•..'here  p^  is  the  pressure  at  infinity,  t  the  time,  J>  the  density  of  fluid. 


a  the  radius  of  the  sphere,  and  ^  the  meridian  angle  measured  from  the 


rear  stagnation  point.  Unfortunately,  however,  the  problem  treated  there 


does  not  exactly  correspond  to  the  present  experiment,  since  the  time- 


isaiVWJiW'JJ' 


dependent  flow  that  can  be  readily  produced  in  the  laboratory  is  a  steady 


uniform  flow  with  sup»*rimposed  traveling  disturbances  with  no  associated 


pressure  fluctuations  in  the  free  stream- 


It  is  not  difficult  to  illustrate  the  possibility  of  such  a  time-dep>endent 


flow.  With  the  Cartesian  coordinates  (X,  Y,  Z),  for  example#  the 


momentum  equations  are  satisfied  by  taking  a  flow  field 


V  =  tj  ,  V  =  f(x-0  t),  V  =  O,  p=  constant,  (2) 
X  o  y  o  z  ' 


which  represent  a  steady  flaw  of  velocity  on  which  a  two-dimensional 


traveling  disturbance  is  superimposed.  No  change  in  pressure  is  produced 


by  the  disturbance.  Since,  however,  the  disturbance  is  rotational  having 


a  component  of  vorticity  in  z-direction,  the  vortex  filaments  are  stretched 


by  the  introduction  of  a  three-dimensional  body  such  as  a  sphere.  This 


means  that  the  "compensating  flow"  introduced  to  satisfy  the  boundary 


conditions  on  the  sphere  must  be  also  be  rotational.  It  does  not  seem  to 


be  impossible  to  find  such  a  compensating  flow,  but  in  view  of  the  limited 


range  of  application  of  the  present  theoretical  calculation,  it  is  doubtful 


whether  the  solution  warrants  the  extra  complication  involved. 


In  order  to  simplify  the  calculations  it  is  assumed  that  the  disturbance 


is  fufficiently  small,  irrotational  and  traveling  v.'ith  a  velocity  equal  to  that 


of  the  steady  uniform  flow.  Using  spherical  polar  coordinates  (r,  X,  > 


such  that 


X  -  r  cos  J?  ,  Y  =  rcin^  cos^,  Z  ~  min  2?  sin  ^  ,  (3) 


the  velocity  potential  is  assumed  in  the  form 


1=  +  5,  + 


(4) 


whei 


=Uo»'(  I  +  2  -^0 


is  the  velocity  potential  for  the  steady  flov/  of  velocity  jjast  a  sphere 

of  radius  a  with  center  at  origin.  The  second  term 

i>,  =  SaU6sinft.(rcoS?9-L/ot;i^r4rs/ni9) 

is  the  velocity  potential  for  the  superimposed  disturbance  traveling  with 

the  phase  velocity  U  and  wave  number  k,  and  finally 

St  Cl^*^ 

=  ea.u.  £  c„.i^  P^ffosi?) 

is  the  velocity  pc  antial  for  the  compensating  flow,  g  is  a  non- 
dimensional  constant,  is  the  modified  Bessel  function  of  the  first  kind, 
P„  the  Legendre  polynomial  of  order  n,  and  C^j  the  non-dimensional 
constants  (n  =  1,  2,  3,  .  .  . )  to  be  determined  by  applying  the  boundary 
conditions  on  the  sphere.  If  6  is  assumed  to  be  small  compared  to  unity, 
the  disturbance  represented  by  gives  a  velocity  fluctuation  of  order 

for  finite  values  of  j8.rsin  2?  ,  and,  when  superimpos  on  the 

steady  uniform  flow,  produces  a  change  in  static  pressure  of  order 
which  is  second  order  and  is  considered  negligible.  Although  the  velocity 
fluctuation  increase  indefinitely  as  i^rs'inJP  tends  to  infinity,  this  does  not 
seem  to  affect  seriously  the  results  provided  that  the  consideration  is 
limited  to  the  region  close  to  the  sphere  and  an  additional  assumption  is 
made  that  ka  is  small  compared  to  unity,  namely,  the  radius  of  the  sphere 
is  small  compared  to  the  wavelength  of  the  disturbance.  Fig.  15  shows 
the  instantaneous  streamlines  (with  equi-distant  values  of  the  Stokes  stream 
function  )  for  the  velocity  potential  for  a  value  ^  = 


Under  the  assumption  that  ka  is  small,  the  normal  velocity 
produced  by  the  disturbance  on  the  sphere  is  given  by  expanding  eq.  (6) 
and  keeping  the  first  two  terms 

(^)  ^  =  C ft<LU,{coS2? COS  kudr  i-ikoLiScDS^  - ijsin 4 Uoij 

which  must  be  cancelled  by  the  normal  velocity  calculated  from  CdJ./dr)^  . 
The  constants  in  the  expression  for  ^2.  then  determined  as 

C,  =  ;^4a_cos6U,'fc,  Ci=54V5m4(jot,  C^-C^- . 

The  pressure  on  the  sphere  is  obtained  from  the  equation 

+  ^  -  fltl 

S  2^  ^  ^  ^  (8) 


^  ”vr  Bfi/rsA.  r  ^  vr 

=  ^U/siij^zP[l  +  le^ika^Ccos  jkoicosfisf/iliU»tj^ 

at  vatA-^'^V.  TW  Jr^A, 

=  —  e  <3L  Uo^f^COi  k  Uot  2?  COSjfisin  kilo  tj* 

Pq  being  the  pressure  at  infinity.  The  pressure  on  the  sphere  is  thus 


given  by 


=|-AU/(lti««a.a»*o.t)  -BeAVU/sin^ U.t  (<,) 


where 


y\=±(-?  1 9cos^(?)  (10) 

4-  /  ^ 

If  the  velocity  V  and  the  acceleration  dV/dt  at  the  location  of  the  sphere 
center  (r  =  o)  are  introduced  by 

y=U.  +-(^)r^  =0.{i+g*«-«s^t/.t) 

^^-aA^AU^shiU^t  J  ,„, 

aZ 
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the  pressure  on  the  sphere  is  expressed  in  the  form 


P-P. 


dt 


(12) 


Eq.  (12)  indicates  -nat  the  pressure  on  the  sphere  consists  of  two  terms. 


the  one  being  proportional  to  the  instantaneous  dynamic  pressure  (1/2)  j^V  , 
and  the  other  proportional  to  the  time  derivative  of  velocity  jDa(dV/dt).  The 
coefficient  A  o^  the  first  term  is  the  same  as  for  the  classical  solution  (1), 
but  the  coefficient  B  of  the  second  term  is  considerably  different. 

On  writing 

y  -  Uo  ^  (13) 


where  &  is  the  meridian  angle  measured  from  the  forward  stagnation 
point,  and  u  is  the  velocity  fluctuation,  (12)  and  (10)  are  written  in  the 
form  j  I 

i-  AUli.  + 


i-9  cos^e)  cO50C-5-f-S’<^i>s*0)\ 


(14) 


5.  DISCUSSION 


In  the  theoretical  calculation  it  was  assumed  that  £  and  ka  are 
small  compared  to  unity,  namely,  the  amplitude  of  velocity  fluctuation  is 
small  compared  to  the  mean  velocity  on  the  jet  axis,  and  the  radius  of 
sphere  is  small  compared  to  the  wavelength  of  velocity  fluctuation.  These 
conditions  were  not  exactly  satisfied  in  the  experiment,  where  the 
maximum  values  of  £  and  ka  are  0.40  and  0.68,  respectively.  Moreover, 
the  observed  phase  velocity  of  velocity  fluctuation  was  only  2/3  of  the  mean 
velocity,  whereas  in  the  theory  the  phase  velocity  is  assumed  to  be  equal 


\ 
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to  the  mean  velocity. 

In  spite  of  these  circumstances,  however,  the  experimental  vaT-ies 
of  the  coefficients  A  and  B  agree  fairly  well  with  theoretical  results 
for  the  sphere  location  at  x  =  6  cm  (Figs.  13  and  14).  The  agreement  is 
net  as  good  for  the  location  at  x  =  3  cm  (Figs.  11  to  12),  and  this  may  be 
due  to  the  undesirable  variation  of  the  velocity  fluctuation  amplitude  around 
that  location. 

As  seen  from  Fig.  2,  the  observed  velocity  fluctuation  is  not  exactly 
simple  sinusoidal,  but  a  slight  distortion  in  the  wave  form  is  apparent  that 
results  in  higher  harmonics.  The  distortion  increases  as  the  fundamental 
frequency  decreases.  This  seems  to  account  for  the  fact  that  the  root- 
mean-square  residual  error  of  representing  the  observed  pressure  fluctua¬ 
tion  by  two  terms  amounts  to  10.  2,  9.7  and  6.  0  per  cer.t  of  the  fluctuation 
amplitude  for  the  fundamental  frequency  f  =  300  Hz  (x  =  3  cm),  350  Hz 
(x  =  6  cm)  and  450  Hz  (x  =  3  and  6  cm),  respectively. 

No  measurement  on  pressure  was  made  on  the  rear  side  of  the  sphere, 
where  the  separation  of  flow  is  expected  to  modify  the  pressure  distribution. 
The  effect  of  separation  can  be  traced  in  the  pressure  distribution  for 
meridian  angle  greater  than  0  =  70®,  where  the  experimental  value  of 
A  begins  to  deviate  from  the  theoretical  curve  (Figs.  13  and  14).  On  the 
other  hand,  the  experimental  value  of  B  agrees  fairly  well  with  the 
theoretical  result  up  to  as  far  as  0  =  90®  (Figs.  13  and  14).  It  is  not 
certain  whether  the  agreement  is  real  or  fortuitous. 
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6-  CONCLUSIONS 


Measurements  of  the  instantaneous  values  of  the  surface  pressure 
were  carried  out  on  a  small  sphere  in  a  periodic  pulsating  flow  and  the 
experimental  values  agreed  moderately  well  with  a  concurrently  developed 
inviscid  theory,  it  is  apparent  from  the  above  theoretical  and  experimental 
work  that  similar  calculations  can  be  performed  on  bodies  other  than 
spheres  or  even  on  a  combination  of  bodies.  By  judicious  choice  of 
such  bodies  points  may  be  found  on  the  surface  where  the  coefficients 
A  =  B  =  O  so  the  measu»-ed  pressure  has  no  contribution  from  the 
acceleration.  Pressure  transducers  at  such  points  or  appropriately 
coupled  to  pressure  holes  at  such  points  would  indicate  static  pressure 
fluctuations  in  the  flow  corresponding  to  the  location  of  the  body  but  with 
the  body  absent  quite  similarlrly  to  rtatic  probes  used  in  steady  flows. 

The  principal  difference  is  that  two  terms  must  be  cancelled  instead  the 
usual  one  in  steady  flow. 
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Fig.  1  Pulsating  flow  generator 
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